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ABSTRACT 

The Milky Way contains several distinct old stellar components that provide a fossil 
record of its formation. We can understand their spatial distribution and kinematics in 
a hierarchical formation scenario by associating the proto-galactic fragments envisaged 
by Searle and Zinn (1978) with the rare peaks able to cool gas in the cold dark matter 
density field collapsing at redshift z > 10. We use hierarchical structure formation sim- 
ulations to explore the kinematics and spatial distribution of these early star-forming 
structures in galaxy haloes today. Most of the proto-galaxies rapidly merge, their stel- 
lar contents and dark matter becoming smoothly distributed and forming the inner 
Galactic halo. The metal-poor globular clusters and old halo stars become tracers of 
this early evolutionary phase, centrally biased and naturally reproducing the observed 
steep fall off with radius. The most outlying peaks fall in late and survive to the present 
day as satellite galaxies. The observed radial velocity dispersion profile and the local 
radial velocity anisotropy of Milky Way halo stars are successfully reproduced in this 
model. If this epoch of structure formation coincides with a suppression of further 
cooling into lower sigma peaks then we can reproduce the rarity, kinematics and spa- 
tial distribution of satellite galaxies as suggested by Bullock et al. (2000). Rcionisation 
at z = 12 ± 2 provides a natural solution to the missing satellites problem. Measuring 
the distribution of globular clusters and halo light on scales from galaxies to clusters 
could be used to constrain global versus local reionisation models. If reionisation oc- 
curs contemporary, our model predicts a constant frequency of blue giobulars relative 
to the host halo mass, except for dwarf galaxies where the average relative frequencies 
become smaller. 

Key words: methods: N-body simulations - methods: numerical - dark matter — 
galaxies: haloes — galaxies: clusters: general, globular clusters 



1 INTRODUCTION 

The Milky Way is a typical bright spiral galaxy. Its disk 
of stars and gas is surrounded by an extended halo of old 
stars, globular star clusters and a few dark matter domi- 
nated old satellite galaxies. For the past 30 years two com- 
peting scenarios for the origin of galaxies and their stellar 
components have driven much observational and theoretical 
research. Eggen, Lynden-Bell and Sandage (1962) proposed 
a monolithic collapse of the Galaxy whilst Searle and Zinn 
(1978) advocated accretion of numerous proto-galactic frag- 
ments. 

Enormous progress has been made in understand- 
ing the structure and origin of the Milky Way, as 
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well as defining a standard cosmological model for 
structure formation that provides us with a framework 
within which to understand our origins IFecblcs 198|; 
iFreeman and Bland-Hawthorn 20021) . Hierarchical growth 
of galaxies is a key expectation within a Universe whose mass 
is dominated by a dark and nearly cold particle (CDM), yet 
evidence for an evolving hierarchy of merging events can 
be hard to find, since much of this activity took place over 
10 billion years ago. The origin of the luminous Galaxy de- 
pends on the complex assembly of its ~ IO^^A/q dark halo 
that extends beyond 200 kpc, and on how stars form within 
the first dark matter structures massive enough to cool gas 
to high densities ( White and Rccs 1978). 

The Galactic halo contains about 100 old metal 
poor globular clusters (i.e. Forbes et al. 2000) each 
containing up to 10® stars. Their spatial distribution 
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falls off as r~'^'^ at large radii and half the globu- 
lars lie within 5 kpc from the centre of the Galaxy 
ijStrader et al. 2004) . There is no evidence for dark mat- 
ter within the globular clusters today jPeebles 1984t 
iMoore 1996h . The old stellar halo population has a similar 
spatial distribution and a total luminos ity of 10** — W^Lq 
jMaiewski et al. 20001: llvezic et al. 2000ll . The stellar popu- 
lations, ages and m etallicities of these compo nents are very 
similar ijFreeman a nd Bland-Ifa wthorn 2002h . 

Also orbiting the Galaxy are several tiny spheroidal 
satellite galaxies, each containing an old population of stars, 
some showing evidence for more recent star-formation in- 
dicating that they can hold on to gas for a Hubble time 
(Gallaaiicr and Wvsc 1994; Grcbcl ct al. 2003). Half of the 
dwarf satellites lie within 85 kpc, have luminosities in the 
range 10^ — 10* L© and are surrounded by dark haloes at 
least 50-200 times as massive as their baryonic components 
iMateo 1998t) . Cold dark matter models have had a no- 
toriously hard time at reconciling the observed low num- 
ber of satellites with the predicted s t eep mass function 
of dark haloes l|Kauffmann et al 19931: iMoore et al. 19991 
iKlvpin et al. 199^ . 

We wish to explore the hypothesis that cold dark mat- 
ter dominates structure formation, the haloes of galaxies 
and clusters are assembled via the hierarchical merging and 
accretion of smaller progenitors (e.g. Lacey and Cole 1993). 
This process violently causes structures to come to a new 
equilibrium by redistributing energy among the collision- 
less mass components. Early stars formed in these progeni- 
tors behave as a collisionless system just like the dark mat- 
ter particles in their host haloes, and they undergo the 
same dynamical processes during subsequent mergers and 
the buildup of larger systems like massive galaxies or clus- 
ters. 

In a recent study, Diemand et al. (2005) used cos- 
mological N-body simulations to explore the distribution 
and kinematics in present-day CDM haloes of dark mat- 
ter particles that originally belonged to rare peaks in the 
matter density field. These properties are particularly rel- 
evant for the baryonic tracers of early CDM structures, 
for example the old stellar halo which may have originated 
from the early disru ption of numerous dwarf proto-galaxies 
jBullock et al. 2000h. the old hal o globular clusters and also 
giant ellipticals jGao et al. 2004ll . 

Since rare, early haloes are strongly biased towards 
overdense regions (e.g. Sheth and Tormen 1999), i.e. towards 
the centers of larger scale fluctuations that have not col- 
lapsed yet, we might expect that the contribution at z = 
from the earliest branches of the merger tree is much more 
centrally concentrated than the overall halo. Indeed, a "non- 
linear" peaks biasing has been discussed by previous author s 
jMoore et al. 19981 : Iwhite and Springel 20001 : iMoore 200ll) . 
Diemand et al. (2005) showed that the present-day distri- 
bution and kinematics of material depends primarily on 
the rareness of the peaks of the primordial density fluc- 
tuation field that the selected matter originally belonged 
to, i.e. when selecting rare density peaks above v<j{M, z) 
[where o{M, z) is the linear theory rms density fluctuations 
smoothed with a top-hat filter of mass M at redshift z\ , their 
properties today depend on u and not on the specific values 
of selection redshift z and minimal mass M. 

In the following section of this paper we discuss a 



model for the combined evolution of the dark and old 
stellar components of the Galaxy within the fra mework 
of the ACDM hierarchical model (|Peebie7T983) . Many 
previous studies have motivated and touched upon as- 
pects of this work but a single formation scenario for 
the above components has not been explored in de- 
tail and compare d w ith data llKauffmann et al. 1995 ; 
Bullock et al. 20001: ICote et al. 20001: iMoore 2001 : 
^ree^a^^^^Kand-Hawtiiprr^OOj : [Benson e^d^002 ; 
Cote et al. 20021: ISp mervillc ct a l. 200S : 

Kravtsov et al. 20041 : iKravtsov and Gnedin 20051 We 
assume proto-galaxies and globular clusters form within 
the first rare peaks above a critical mass threshold that 
can allow gas to cool and form stars in significant numbers 
(typically at z « 12). We assume that shortly after the 
formation of these first systems, the universe reionises, 
perhaps by these first proto-galaxies, suppressing further 
formation of cosmic structure until later epochs. We use the 
N-body simulations to trace the rare peaks to z = 0. Most 
of these proto-galaxies and their globular clusters merge 
together to create the central galactic region. In Section 
3 we will compare the spatial distribution and orbital 
kinematics of these tracer particles with the Galactic halo 
light and old metal poor globular clusters. We will see 
that a small number of these first stellar systems survive 
as dark matter dominated galaxies. We will compare their 
properties with the old satellites of the Galaxy in Section 4. 



2 THE FIRST STELLAR SYSTEMS 

We propose that 'ordinary' Population II stars and glob- 
ular clusters first appeared in significant numbers at red- 
shift > 12, as the gas within protogalactic haloes with virial 
temperatures above lO'^K (corresponding to masses com- 
parable to those of present-day dwarf spheroidals) cooled 
rapidly due to atomic processes and fragmented. It is this 
'second generation' of subgalactic stellar systems, aided per- 
haps by an earlier generation of metal-free (Population 
HI) stars and by their remnant black holes, which gen- 
erated enough ultraviolet radiation to reheat and reion- 
ize most of the hydrogen in the Universe by a redshift 
z = 12, thus preventing further accretion of gas into the 
shallow potential wells that collapsed later. The impact of 
a high redshift UV background o n structure formation has 
been invoked by se v eral authors fHaardt and Madau 199J ; 
BuUock et al. 20001: iMoore 200ll: .Barkana and Loeb 20ol : 
TuUv et al. 20021: iBenson et al. 2002^ to explain the flatten- 
ing of the faint end of the luminosity function and the miss- 
ing satellites problem within our Local Group. Here we use 
high resolution numerical simulations that follow the full 
non-linear hierarchical growth of galaxy mass haloes to ex- 
plore the consequences and predictions of this scenario. 

Dark matter structures will collapse at different times, 
depending on their mass, but also on the underlying larger 
scale fluctuations. At any epoch, the distribution of masses 
of collapsed haloes is a steep power law towards low masses 
with n(m) oc m^^ . To make quantitative predictions we 
calculate the non-linear evolution of the matter distribu- 
tion within a large region of a ACDM Universe. The en- 
tire well resolved region is about 10 comoving megaparsecs 
across and contains 61 million dark matter particles of mass 
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Figure 1. The high redshift and present day mass distribution 
in a region that forms a single galaxy in a hierarchical cold dark 
matter Universe. The upper panel shows the density distribution 
at a redshift z = 12 from a region that will form a single galaxy 
at z = (lower panel). The blue-pink colour scale shows the 
density of dark matter whilst the green regions show the particles 
from proto-galaxies with virial temperature above 10'' K that 
have collapsed at this epoch. These peaks have masses in the range 
10® — 10^" Mq. The lower panel shows same mass distribution at 
2 = 0. Most of the rare peaks are located towards the centre of 
the galaxy today. The squares in both panels indicate those first 
objects that survive the merging process and can be associated 
with the visible satellite galaxies today orbiting within the final 
galactic mass halo. Most of the subhaloes stay dark since they 
collapse later after reionisation has increased the Jeans mass. 



5.85 X 10^ and force resolution of 0.27 kpc. This region 
is embedded within a larger 90 Mpc cube that is simulated 
at lower resolution such that the large scale tidal field is 
represented. Figure 1 shows the high-redshift and present- 
day mass distribution of a single galaxy mass halo taken 
from this large volume. The rare peaks collapsing at high 



redshift that have had sufficient time to cool gas and form 
stars, can be identified, followed and traced to the present 
day. Because small fiuctuations are embedded within a glob- 
ally larger perturbation, the small rarer peaks that collapse 
first are closer to the centre of the final potential and they 
preserve their locality in the present day galaxy. The strong 
correlation between initial and final position results in a sys- 
tem where the oldest and rarest peaks are spatially more 
concentrated than less rare peaks. The present day spatial 
clustering of the material that was in collapsed structures at 
a higher redshift only depends on the rarity of these peaks 
jPiemand et a,]. 200,^^. 

Our simulation contains several well resolved galactic 
mass haloes which we use to trace the evolution of progeni- 
tor haloes that collapse at different epochs. The first metal 
free Population III stars form within minihaloes already col- 
lapsed hy z > 25, where gas can cool via roto-vibrational 
levels of H2 and contract. Their evolution is rapid and lo- 
cal metal enrichment occurs from stellar evolution. Metal- 
poor Population II stars form in large numbers in haloes 
above Mh ^ 10** [(1 -f z)/10]"^/^ Mq (virial temperature 
10" K), where gas can cool efficiently and fragment via ex- 
citation of hydrogen Lya. At z > 12, these correspond to 
> 2.5 a peaks of the initial Gaussian overdensity field: most 
of this material ends up within the inner few kpc of the 
Galaxy. Within the « IMpc turn-around region, a few hun- 
dred such protogalaxies are assembling their stellar systems 
l|Kravtsov and Gnedin 20 05^. Typically 95% of these first 
structures merge together within a timescale of a few Gyrs, 
creating the inner Galactic dark halo and its associated old 
stellar population. 

With an efficiency of turning baryons into stars and 
globular clusters of the order /* = 10% we successfully re- 
produce the total luminosity of the old halo population and 
the old dwarf spheroidal satellites. The fraction of baryons 
in dark matter haloes above the atomic cooling mass at 
redshift 12 exceeds fc — 1%. A normal stellar population 
with a Salpeter-type initial mass function emits about 4,000 
hydrogen-ionizing photons per stellar baryon. A star for- 
mation efficiency of 10% therefore implies the emission of 
4, 000 X ft X fc a, few Lyman-continuum photons per 
baryon in the Universe. This may be enough to photoionize 
and drive to a higher adiabat vast portions of the inter- 
galactic medium, thereby quenching gas accretion and star 
formation in nearby low-mass haloes. 



3 CONNECTION TO GLOBULAR CLUSTERS 
AND HALO STARS 

The globular clusters that were once within the merging 
proto-galaxies are so dense that they survive intact and will 
orbit freely within the Galaxy. The surviving proto-galeixies 
may be the precursors of the old satellite galaxies, some 
of which host old globular clusters such as Fornax, whose 
morphology and stellar populations are determined by on- 
going gravitational and hydrodynamical interactions with 
the Milky Way (e.g. Mayer et al. 2005). 

Recent papers have attempted to address the origin of 
the spatial distribution of globular clusters (e.g. Parmentier 
and Grebel 2005, Parmentier et al. 2005). Most compelling 
for this model and one of the key results in this paper, is 
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Figure 2. The radial distribution of old stellar systems com- 
pared with rare peaks within a z = LCDM galaxy. The thick 
blue curve is the total mass distribution today. The labeled green 
curves show the present day distribution of material that col- 
lapsed into 1,2,2.5,3 and 3.5a peaks at a redshift z = 12. The 
circles show the observed spatial distribution of the Milky Way's 
old metal poor globular cluster system. The dashed line indicates 
a power law p(r) oc r~^'^ which represents the old halo stellar 
population. The squares show the radial distribution of surviv- 
ing 2.5(7 peaks which are slightly more extended than the overall 
NFW like mass distribution, in good agreement with the observed 
spatial distribution of the Milky Way's satellites. 

that we naturally reproduce the spatial clustering of each of 
these old components of the galaxy. The radial distribution 
of material that formed from > 2.5a peaks at 2 > 12 now 
falls off as p{r) oc r~^'^ within the Galactic halo - just as 
the observed old halo stars and metal poor globular clusters 
(cf Figure 2). Cosmological hydrodynamical simulations are 
also begining to attain the resolution to resolve the forma- 
tion of the old stellar haloes of galaxies ( Abadi et al. 2005) . 
Because of the steep fall off with radius, we note that we do 
not expect to find any isolated globular clusters beyond the 
virial radius of a galaxy ^. 

These first collapsing structures infall radially along fil- 
aments and end up significantly more flattened than the 
mean mass distribution. They also have colder velocity dis- 
tributions and their orbits are isotropic in the inner halo and 
increasingly radially anisotropic in the outer part. Material 
from these rare peaks has f3 = 1 — {v^/v^) « 0.45 at our 
position in the Milky Way, in remarkable agreement with 
the recently measured anisotropy and velocity dispersion 
of halo stars (Chiba and Beers 2000; Battaglia ct al. 2005; 
iThom et ai~2005ill . Diemand et al. (2005) show that the ra- 

^ The probability to have one isolated old globular cluster outside 
of the virial radius of a Milky Way like galaxy is only 3% in our 
model. 



dial distribution of rarer peaks is even more highly biased - 
thus the oldest population III stars and their remnant black 
holes are found mainly within the inner kpc of the Galaxy, 
falling off with radius steeper than r~*. 

The observational evidence for tidal stellar streams from 
globular clusters suggests that they are n ot embedded within 
extended dark matter structures today jMoore 19961) . This 
does not preclude the possibility that the globular clusters 
formed deep within the central region of 10* M0 dark haloes 
which have since merged together. (Massive substructure 
within the inner ~ 2Q%R^iriai of galactic mass haloes is 
tidally disrupted i.e. Gihgna et al. 1998.) This is what we 
expect within our model which would leave the observed 
globulars freely orbiting without any trace of the original 
dark matter component. However, it is possible that the 
most distance halo globulars may still reside within their 
original dark matter halo. If the globular cluster is located 
at the center of the CDM cusp, then observations of their 
stellar kinematics may reveal rising dispersion profiles. If the 
globular cluster is orbiting within a CDM mini-halo then 
we would expect to see symmetric tidal streams orbiting 
within the potential of the CDM substructure halo rather 
than being stripped by the Galaxy. 



4 CONNECTION TO SATELLITE GALAXIES 
AND THE MISSING SATELLITES 
PROBLEM 

The remaining ~ 5% of the proto-galaxies form sufficiently 
far away from the mayhem that they fall into the assem- 
bling galaxy late (2 ~ 1 — 2, about one Gyr after the for- 
mation of the inner Galaxy at z ~ 5). This leaves time 
to enha nce their g/Fe element r a tios from Type II super- 
novae llWvse and Gilmore 1988t IWvse and Gilmore 1991: 
IPritzl et al. 200R^ . Recent studies including chemical mod- 
eling of this process support this scenario (e.g. Robertson et 
al. 2005, Font et al. 2005). 

The proto-galaxies highlighted with boxes in Figure 1 
are those few systems that survive until the present epoch - 
they all form on the outskirts of the collapsing region, end- 
ing up tracing the total mass distribution as is also observed 
within the Milky Way's and M31's satellite systems. Each of 
our four high resolution galaxies contains about ten of these 
surviving proto-galaxies which have a radial distribution 
that is slightly shallower than that of the total mass distribu- 
tion but more concentrated than the distribution of all sur- 
viving (or 2 = mass selected) subhalos (Figures |21 and . 
This is consistent with the spatial distribution of surviving 
satellites i n the Milky Way and in other nearby galaxies in 
the 2dF (van den Bosch et al. 20051: ISales fc Lambas 20051) 
and DEEP2 samples JCoil et al. 20051) and with g alaxy 
groups like NGC5044 llFaltenbacher fc Mathews 20051) . 

Figure 3 shows the distribution of circular velocities of 
the Local Group satellites compared with these rare proto- 
galaxies that survive until the present day. The Local Group 
circular velocity data are the new data from Maccio' et al. 
(2005b) where velocity dispersions have been converted to 
peak circular velocities using the results of Kazantzidis et 
al. (2004). The total number of dark matter substructures 
is over an order of magnitude larger than the observations. 
Reionisation and photo-evaporation must play a crucial role 
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Figure 3. The cumulative velocity distribution function of ob- 
served Local Group satellites and predicted dark matter substruc- 
tures. The red squares show the observed distribution of circu- 
lar velocities for the Local Group satellites. Cold dark matter 
models predict over an order of magnitude more dark matter 
substructures than are observed within the Galactic halo (blue 
curve). The black solid curve show the cumulative velocity dis- 
tribution of present day surviving substructure haloes that were 
the rare > 2.5(t peaks identified ai z = 12, before they entered a 
larger mass halo. The dashed curve shows the same objects but at 
2 = 1.1 before they entered a larger system in the hierarchy. Their 
similarity shows that little dynamical evolution has occurred for 
these objects. 



in suppressing star formation in less rare peaks, thus keeping 
most of the low mass haloes that collapse later devoid of 
baryons. The surviving population of rare peaks had slightly 
higher circular velocities just before accretion (at z ~ 1, 
dashed line in Figure 3 - see Kravtsov et al. 2004), tidal 
stripping inside the Galaxy halo then reduced their masses 
and circular velocities and they match the observations at 
z = Q. Dissipation and tidally induced bar formation could 
enable satellites to survive even closer to the Galactic centre 
(Maccio' et al. 2005a). 

Likewise to the radial distribution, the kinematics of 
the surviving visible satellite galaxies resembles closely the 
one of the dark matter while the same properties for all 
the surviving subhalos differ (Figures |^ and [SJ . Within the 
four high resolution CDM galaxy haloes our 42 satellite 
galaxies have average tangential and radial velocity dis- 
persions of 0.70±0.08Vc,max and 0.56±0.07Vc,max respec- 
tively, i.e. P = 0.26 ± 0.15 (the errors are one sigma Pois- 
sion uncertainties). These values are consistent with those 
of the dark matter particles: crtan ~ 0.66Vc,max, cTrad ~ 
0.55Vc,max and 13 = 0.30; the hint of slightly larger disper- 
sions of the satellites are consistent with their somewhat 
larger radial extent. In the inner part our model satellite 
galaxies are hotter than the dark matter background, es- 
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Table 1. Present-day properties of the four simulated galaxy 
haloes. The columns give halo name, virial mass, virial radius, 
peak circular velocity, and radius to the peak of the circular 
velocity curve. The virial radius is defined to enclose a mean 
density of 98.6 times the critical density. The mass resolution 
is 5.85 X 10^ Mq and the force resolution (spline softening length) 
is 0.27 kpc. For comparison with the Milk Way these halos were 
rescaled to a peak circular velocity of 195 km/s. In SPH simu- 
lations of the same halos we found that this rescaling leads to a 
local rotation speed of 220 km/s after the baryonic contraction 
iMaccio' et al. 2005^ . The rescaled virial radii and virial masses 
are given in the last two columns. 
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pecially in the tangential component: Within 0.3 

find CTrad.GALs/o-rad.DM = 0.69 Vcmax /0.62 Vc.max = l.H and 
ITtan.GALs/o-tan.DM = 0.95 Vc, max /0. 76 14, max = 1-25. This is 

consistent with the observed radial velocities of Milky Way 
satellites. For the inner satellites also the tangential mo- 
tions are know (with large uncertainties however) (e.g. Ma- 
teo 1998; Wilkinson & Evans 1999) and just as in our simple 
model they are larger than the typical tangential velocities 
of dark matter particles in the inner halo. 

The total (mostly dark) surviving subhalo popula- 
tion is more extended and hotter than the dark mat- 
ter while the distri bution of orbits (i.e. 13) is similar 
JPiemand et al. 20041) . For the 2237 subhalos within the four 

galaxy haloes find atan = 0.84Vc,max, crtan = 0.67yc,max 

and (3 = 0.21, i.e. there is a similar velocity bias rel- 
ative to the dark matter in both the radial and tan- 
gential components and therefore a similar anisotropy. In 
the inner halo the differences between dark matter parti- 
cles and subhaloes are most obvious: Within 0.3 r-vir we 

find arad,SUB/o-rad,DM = 0.91K,max/0.62 Vc.max = 1.47 and 
0-tan,SUB/o-tan,DM = 1 .21 Vc,max/0.76Vc,max = 1-59. SubhaloS 

tend to avoid the inner halo and those who lie near the cen- 
ter at 2 = move faster (both in the tangential and radial 
directions) than the dark matter particles, i.e. these inner 
subhalos have large orbital energies and spend most of their 
time further away from the center (Figures |1] |^ and |S| see 
also Diemand et al. 2004). 



5 SUMMARY 

We have a implemented a simple prescription for proto- 
galaxy and globular cluster formation on to a dissipationless 
CDM N-body simulation. This allows us to trace the kine- 
matics and spatial distribution of these first stellar systems 
to the final virialised dark matter halo. We can reproduce the 
basic properties of the Galactic metal poor globular cluster 
system, old satellite galaxies and Galactic halo light. 

The spatial distribution of material within a virialised 
dark matter structure depends on the rarity of the peak 
within which the material collapses. This implies a degener- 
acy between collapse redshift and peak height. For example, 
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Figure 5. Radial velocities of satellite galaxies (filled squares), all (dark and luminous) subhalos (open squares) and dark matter particles 
(small points) in four CDM galaxy halos. The solid lines are the radial velocity dispersion of the dark matter plotted with positive and 
negative sign. All quantities are in units of the virial radius and maximum of the circular velocity of the host halos. For comparison with 
the observed satellite galaxies around the Milky Way (filled circles) from (Matoo 1998; Wilkinson & Evans 1999) the simulated halos 
were rescaled (see Table l4l . The observed and modeled satellite galaxies have similar radial velocities as the dark matter particles while 
those of the dark subhalos are larger, especially in the inner part. 



3 sigma peaks collapsing at redshift 18 and 10 will have the 
same final spatial distribution within the Galaxy. However 
this degeneracy can be broken since the mass and number 
of peaks are very different at each redshift. In this example 
at redshift 18 a galaxy mass perturbation has 700 collapsed 
3 sigma halos of mass 6 x 10^ Mq, compared to 8 peaks of 
mass 4 x W^Mq 

The best match to the spatial distribution of globular 



clusters and stars comes from material that formed within 
peaks above 2.5 a. We can then constrain the minimum 
mass/redshift pair by requiring to match the observed num- 
ber of satellite galaxies in the Local Group (Figure |7|l. If 
protogalaxies form in early, low mass 2.5 a peaks the result- 
ing number of luminous satellites is larger as when they 
form later in heavier 2.5 a peaks. We find that efficient 
star formation in halos above about 10* up to a red- 
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Figure 6. Tangential velocities of satellite galaxies (filled squares), all (dark and luminous) sublialos (open squares) and dark matter 
particles (small points) in four CDM galaxy halos. The lines are the tangent ial velocity dispersion of the dark matter (solid) and the 
circular velocity (dashed). The four satellite galaxies iWilkinson fc: E vans 19931) give tangential velocities (i.e. from inside out: LMC/SMC, 
Ursa Minor, Sculptor and Draco) are plotted with fi lled circle s. The open triang les with error bars show HST pr oper motion data for 
(from the inside out) Ursa Minor JPiatek et al. 200fi) . Carina iPiatek et al. 2n03h and Fornax iPiatek et al. 20o3) . The units are as in 
Figure |5] The observed and modeled inner satellite galaxies (and also the dark inner subhalos) have larger typical tangential velocities 
than the dark matter particles in the same regions. 



shift z = ll.Sl^'g matches these constraints. The scatter in 
redshift is due to the different best fit redshifts found in our 
individual galaxy haloes. After this epoch star formation 
should be suppressed in small halos otherwise a too large 
number of satellites and a too massive and too extended 
spheroid of population II stars is produced. The minimum 
halo mass to form a protogalaxy inferred from these two 



constraints corresponds to a minimal halo virial tempera- 
ture of lO'^K (Figure |7|l, i.e. just the temperature needed 
for efficient atomic cooling. 

This model is general for galaxy formation, but individ- 
ual formation histories may reveal more complexity. Soon af- 
ter reionisation, infalling gas into the already massive galac- 
tic mass halo leads to the formation of the disk and the metal 
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r [r - - J r [r J expansion factor a = (z+l) ' 



Figure 4. Enclosed number of satellite galaxies (solid lines), all 
(dark and luminous) subhalos (dotted) and dark matter particles 
(dashed) in four CDM galaxy halos. The numbers are relative to 
the value within the virial radius. The subhalos are only plot- 
ted out to the viral radius. For comparison with the observed 
satellite galaxies around the Milky Way (dash-dotted lines) from 
iMat eo 1998; Wilkinson & Evans 1999) the simulated halos were 
rescaled (see Table The satellite galaxy distribution is more 
concentrated than the one of the total surviving subhalo popu- 
lation but usually more extended than the dark matter particle 
distribution but there are large differences from one halo to the 
other. Well beyond the virial radius, the numbers of field dwarf 
galaxies that will host stars falls below the mean dark matter 
density. 



enriched population of globular clusters. The first and sec- 
ond generation of stars forming in proto-clusters of galaxies 
will have a similar formation path, but occurring on a more 
rapid timescale. 

We find that the mass fraction in peaks of a given 
a is independent of the final halo mass, except that it 
rapidly goes to zero as the host halos become too small 
to have suf ficiently high a proge nitors (see Figure |H| and 
Table 4 in jPiemand et al. 2005^ ). Therefore, if reiomsa- 
tion is globally coeval throughout the Universe the abun- 
dance of globulars normalised to the halo mass will be 
roughly constant in galaxies, groups and clusters. Further- 
more, the radial distribution of globular clusters relative 
to the host halo scale radius will the same (see Diemand 
et al. 2005). If ra rer peaks reionise galaxy clusters earlier 
JTuUv et aL 2002i) then their final distribution of blue glob- 
ulars will fall off more steeply (relative to the scale radius 
of the host halo) and th ey will be less abundant per virial 
mass jPiemand et al. 20 05). Observations suggest that the 
numbers of old globular clusters are correlated wit h the lu- 
minosity of t he host galaxy ( . McL au ghhn 1999 ; Harris 199lt 
iHarris 20051 : iRhode et al. 20051) . Wide field surveys of the 
spatial distribution of globulars in groups and clusters may 



Figure 7. Minimal halo mass for protogalaxy formation vs. time. 
The dashed line connects minimal mass/redshift pairs which pro- 
duce the number of satellite galaxies observed in the Local Group. 
If efficient protogalaxy formation is allowed below this line our 
model would produce too many satellites. The best match to the 
spatial distribution of globular clusters and stars comes from ma- 
terial that formed within peaks above 2.5 a (dash-dotted). The 
circle with error bars indicates the latest, lowest mass halo which 
is still allowed to form a protogalaxy. The uncertainty in redshift 
z = 11.5^^ 5 is due to scatter in best fit redshift when matching 
the the spatial distribution of globular clusters to our individ- 
ual galaxy halo models at a fixed minimum mass of 10^ Mq. The 
range in minimum masses produces N^at ± V^aat ~ 11 it 3 lumi- 
nous subhalos around an average galaxy. The dotted line shows 
halos with the atomic cooling virial temperature of lO'^K. Our 
inferred minimal mass for efficient protogalaxy formation follows 
the dotted line until z = ll.sjg'^ and rises steeply (like the 2.5 
a line or steeper) after this epoch. 

reveal the details of how and when reionisation occurred 
jForbes et al. 19971; IPu^ia et al. 2001) . 
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